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Abstract 

In  this  paper  we  present  three  physiologically  based  pharmacokinetic  (PBPK)  models  for  the  systemic 
transport  of  trichloroethylene  (TCE),  with  a  focus  on  the  adipose,  or  fat  tissue.  TCE  is  a  widespread 
environmental  contaminant,  and  has  been  shown  to  produce  toxic  effects  in  both  animals  and  humans. 
A  key  characteristic  of  TCE  is  its  tendency  to  accumulate  in  fat  tissue,  which  has  a  major  impact  on  the 
overall  systemic  disposition  of  TCE. 

Here  we  use  PBPK  models  to  predict  the  kinetics  of  TCE  through  the  various  tissues  and  organs, 
including  the  adipose  tissue.  The  first  model  utilizes  the  standard  “perfusion-limited”  compartmental 
model  for  the  fat  tissue,  while  the  second  model  uses  a  “diffusion-limited”  model  to  describe  the  transport 
through  the  adipose  tissue.  Both  of  these  ODE  models  are  based  on  ’’well-mixed”  and  rapid  equilibrium 
assumptions,  and  do  not  take  into  account  the  specific  and  largely  heterogeneous  physiology  of  adipose 
tissue. 

The  third  model  we  discuss  is  a  PBPK  hybrid  model  with  an  axial-dispersion  type  model  for  the  adi¬ 
pose  tissue.  This  PDE-based  model  is  designed  to  capture  key  physiological  heterogeneities  of  fat  tissue, 
including  widely  varying  fat  cell  sizes,  lipid  distribution,  and  blood  flow  properties.  Model  simulations 
demonstrate  that  this  model  may  be  well-suited  to  predict  the  experimental  behavior  of  TCE  in  adipose 
tissue  using  parameter  estimation  techniques. 


1  Introduction 


Many  environmental  and  occupational  toxins  are  lipophilic  in  nature,  and  tend  to  accumulate  in  the  adipose, 
or  fat  tissue  of  animals  and  humans.  These  compounds  vary  widely  in  their  means  of  uptake  into,  disposition 
through,  and  elimination  from  the  body,  as  well  as  in  their  toxic  mechanisms.  A  key  component  in  under¬ 
standing  their  specific  toxicities  is  the  characterization  of  their  kinetics.  Since  adipose  is  a  storage  tissue  for 
these  xenobiotics,  it  is  important  to  ascertain  their  absorption,  distribution  and  accumulation  properties  for 
adipose  tissue. 

Numerous  studies  have  shown  that  a  large  class  of  lipophilic  compounds  accumulates  in  adipose  tissue. 
This  class  includes  organochloride  pesticides  such  as  aldrin,  dieldrin,  endrin,  chlordane,  heptachlor  and 
DDT  [2],  as  well  as  polychlorinated  biphenyls  [2], [36],  polychlorinated  dibenzo-p-dioxins  and  polychlorinated 
dibenzo-furans  [15].  Moreover,  a  variety  of  other  halogenated  organic  compounds  are  known  to  accumulate  in 
fat  tissue  [2].  This  tissue  accumulation  is  related  to  the  lipophilicity  of  the  compounds  and  the  low  perfusion 
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rate  of  adipose  tissue,  although  studies  in  [21]  and  [39]  suggest  that  the  tissue  distribution  of  lipophilic 
compounds  is  also  dependent  on  binding  competition  between  fat  tissue  and  well-perfused  lean  tissues. 

As  a  motivating  example,  we  consider  the  solvent  trichloroethylene  (TCE),  which  is  widely  used  in  in¬ 
dustry  as  a  metal  degreasing  agent.  Humans  come  into  contact  with  TCE  most  often  by  inhaling  TCE 
vapor  or  by  drinking  contaminated  water.  TCE  is  a  highly  lipophilic  volatile  organic  compound,  and  is 
quickly  distributed  into  the  bloodstream.  It  is  partially  eliminated  from  the  body  through  exhaled  air 
and  is  rapidly  metabolized  in  the  liver.  There  are  two  major  metabolic  pathways  for  TCE:  Cytochrome 
P450-mediated  oxidation  and  glutathione  conjugation  [33].  Major  metabolites  include  chloral  hydrate, 
trichloroacetic  acid  (TCA),  dichloroacetic  acid  (DCA),  free  trichloroethanol,  trichloroethanol-glucuronide 
and  S-(l,2-dichlorovinyl)-L-cysteine  (DCVC)  [20], [33],  and  most  are  eliminated  in  the  urine  [50]. 

Unmetabolized  TCE  accumulates  in  adipose  tissue  [50],  but  over  long  periods  of  time  is  not  very  persistent 
in  the  body  because  of  its  rapid  metabolism.  In  a  study  by  Bergman  [8] ,  mice  were  exposed  to  a  dose  of  280 
mg/kg  radiolabeled  TCE  for  ten  minutes,  and  the  tissue  distribution  of  TCE  was  measured  using  whole- 
body  autoradiography.  Within  30  minutes,  there  was  a  large  accumulation  of  TCE  in  the  fat  tissue,  with 
significant  levels  still  present  after  four  hours  and  traces  still  visible  after  eight  hours. 

TCE  and  several  of  its  metabolites  are  known  to  produce  toxic  effects  in  humans.  Acute  exposure 
to  TCE  impacts  the  central  nervous  system,  causing  symptoms  such  as  fatigue,  headaches,  dizziness  and 
drowsiness  [28].  Several  of  its  metabolites  are  suspected  of  causing  toxicity  in  animals  and/or  humans, 
including  TCA,  DCA,  DCVC  and  chloral  hydrate  [13], [40].  Evidence  suggests  that  both  TCA  and  DCA  are 
hepatocarcinogens  [13],  [14]  while  chloral  hydrate  is  associated  with  lung  tumors  in  mice  [29].  Moreover,  a 
recent  study  shows  that  TCE  may  play  a  role  in  the  development  of  kidney  cancer  in  humans  [12],  and  the 
metabolite  DCVC  has  been  linked  to  kidney  tumors  in  rats  [34].  TCE  and  some  of  its  metabolites  also  have 
been  associated  with  developmental  toxicity  in  animals  [43],  including  eye  and  cardiac  malformations  [7]. 

These  lipophilic  environmental  toxins  such  as  TCE  can  have  both  short-term  and  long-term  adverse 
health  effects  to  our  population.  It  is  important  to  assess  the  health  risks  these  chemicals  pose,  so  that  the 
public  can  be  adequately  protected  from  dangerous  compounds  and  the  safer  compounds  can  be  used  with 
greater  confidence.  In  order  to  better  determine  how  xenobiotics  such  as  TCE  cause  toxic  effects  in  the  body, 
the  toxicokinetics  of  these  compounds  must  be  understood.  Physiologically  based  pharmacokinetic  (PBPK) 
models  are  developed  to  describe  the  disposition  of  a  given  chemical  throughout  the  body.  By  characterizing 
the  kinetics  of  the  compound,  we  can  predict  its  transport  inside  the  body  and  work  towards  gaining  a  better 
understanding  of  its  toxic  effects  and  mechanisms. 

The  standard  PBPK  models,  which  make  use  of  “perfusion-limited”  and  “diffusion-limited”  compart¬ 
ments,  are  designed  around  the  assumption  of  rapid  equilibrium  and  uniformity  within  each  compartment  or 
subcompartment.  This  is  often  a  reasonable  assumption  in  physiological  modeling,  but  the  nonhomogeneous 
composition,  structure  and  circulation  properties  of  adipose  tissue  suggest  that  this  type  of  “well-stirred” 
model  may  not  be  sufficient.  An  axial  dispersion  model,  such  as  the  one  derived  by  Roberts  and  Rowland 
[41]  for  hepatic  uptake  and  elimination,  is  designed  to  account  for  such  physiological  heterogeneities. 

In  this  paper  we  develop  toxicokinetic  models  to  describe  the  behavior  of  inhaled  trichloroethylene  as  it 
enters,  flows  through,  and  possibly  accumulates  in  adipose  tissue.  We  formulate  several  physiologically  based 
models,  each  with  a  different  submodel  for  the  adipose  tissue  component.  First  we  derive  the  standard  PBPK 
model  with  perfusion-limited  compartments.  We  then  modify  this  model  by  replacing  the  perfusion-limited 
adipose  tissue  compartment  with  a  diffusion-limited  compartment.  Finally,  we  adapt  the  dispersion  model 
of  Roberts  and  Rowland  to  describe  the  transport  of  TCE  through  adipose  tissue.  This  compartment  is  then 
coupled  with  the  remaining  perfusion-limited  compartments  to  create  a  whole-body  PBPK  model. 
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Standard  PBPK  models 


2.1  A  perfusion-limited  PBPK  model  for  TCE 


Many  PBPK  models  for  lipophilic  compounds  include  compartments  for  tissues  such  as  the  liver,  lungs, 
adipose  tissue,  richly  perfused  and  poorly  perfused  tissues.  These  compartments  often  assume  a  perfusion- 
limited  model,  or  equivalently,  a  flow-limited  model  of  disposition,  meaning  that  the  rate  of  uptake  of  the 
compound  into  the  tissue  is  limited  by  the  blood  flow  rate  to  the  tissue  rather  than  the  rate  of  diffusion 
across  the  cell  membranes  [38].  In  this  case,  the  blood  and  tissue  are  in  equilibrium  and  the  equation  for 
motion  of  a  solute  through  the  constant-volume  tissue  compartment  is  an  ordinary  differential  equation  of 
the  form 


dC 

v—  =  Q{Cin  ~  Cout) ’ 


where  V  is  the  volume  of  the  tissue  in  liters,  C  is  the  concentration  of  compound  inside  the  tissue  (in 
mg/liter),  Q  is  the  blood  flow  rate  to  the  tissue  (in  liters/hour),  and  Cj„  and  Cout  are  the  compound 
concentrations  entering  and  exiting  the  tissue,  respectively.  Examples  of  PBPK  models  for  TCE  and  its 
metabolites  that  use  perfusion-limited  compartments  can  be  found  in  [1],  [24],  [25],  [30],  [48]. 


Here  we  develop  a  standard  PBPK  model  [38]  for  inhaled  TCE  with  flow-limited  compartments  for  the 
kidney,  muscle  tissue,  adipose  tissue,  brain,  liver,  venous  blood,  and  remaining  non-fat  tissue  (see  Figure  1). 
We  assume  uptake  via  inhalation,  with  a  lung  compartment  subdivided  into  the  alveolar  space  and  lung 
blood  subcompartments.  TCE  is  metabolized  in  the  liver,  which  is  modeled  with  Michaelis-Menten  kinetics. 


In  the  lung,  ventilation  is  assumed  to  be  continuous  with  rate  Qp,  and  the  vapor  in  the  alveolar  space  is 
in  rapid  equilibrium  with  the  arterial  lung  blood,  which  has  cardiac  output  rate  Qc.  Moreover,  the  blood/air 
partition  coefficient  is  denoted  by  P&. 


The  variables  used  in  the  lung  compartment  include: 


Cc 

Cx 

Ca 

Cv 

A-i 

Ax 

Al 


Concentration  of  xenobiotic  in  surrounding  air 
Concentration  of  xenobiotic  in  alveolar  space 
Concentration  of  xenobiotic  in  arterial  blood 
Concentration  of  xenobiotic  in  venous  blood 
Amount  of  xenobiotic  inhaled 
Amount  of  xenobiotic  exhaled 
Amount  of  xenobiotic  in  lung. 


In  this  case, 
arterial  blood: 


the  concentration  Cx  in  the  alveolar  air  is  related  linearly  to  the  concentration  Ca  in  the 


Cx 


Ca 

Pb' 


The  rate  of  inhalation  of  xenobiotic  is  given  by  QpCc,  while  the  rate  of  exhalation  is  given  by  QPCX. 
Therefore  we  have  the  following  equations: 


dAj 

dt 

dAx 

dt 

dA^ 

dt 


QPCc 

QPCX  =QP %■ 


QP(Cc-Cx)  +  Qc(Cv  -Ca). 
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Moreover,  the  assumptions  of  the  model  [38]  imply  that  =  0,  so  by  substituting  Cx  =  Ca/Pb  we  obtain 


C0  = 


QpCc  “b  QcCv 

Qc  +  T\ 


Combining  the  perfusion-limited  compartments  with  the  lung  compartment,  we  obtain 

=  Qf(Ca-Cvf ) 


Vv 


dt 

dCv 

dt 

Ca 


v„ 


dCn 


Vt 

vhr 


dt 

dC± 

dt 

dCbr 


dt 


QmCvm  “b  QtCvt  “b  QfCvf  "b  QbrCvbr  ~b  QlCvl  ~b  QkCvk  QcCv 
Q  cP'V  +  QpCc 
Qc  +  t% 

Qm(Ca  ~  Cvm) 

Qt(Ca  ~  Cvt) 

Qbr^C'a  C'vbr') 


rlA 

'VmaxC'vl 

dt 

kM  +  Cvl 

V^ 

1  dt 

=  Ql{Ca  —  Cvl)  ■ 

-S£  1 

e>  N 

-e  1 
£ 

—  Qk(Ca  Cvk ) 

dAi 

dt 

=  QpCc 

dAx 

dt 

—  QpCx- 

V-max  cvl 

kjvi  +  Cvi 


The  subscripts  denote  the  following  specific  tissues: 


V 

Venous  blood 

k 

Kidney 

m 

Muscle 

f 

<£> 

Fat 

br 

<£> 

Brain 

l 

<£> 

Liver 

t 

Remaining  non-fat  tissue. 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 
(9) 

(10) 

(11) 


The  amount  of  solute  metabolized  in  the  liver  is  denoted  by  Aam ,  and  has  units  in  milligrams.  Con¬ 
stants  in  the  model  include  the  Michaelis-Menten  constant  kM  (mg/liters)  and  the  metabolic  constant  vmax 
(mg/  (hour)). 

Volumes  (in  liters)  of  specific  tissues  are  denoted  by  V,  concentrations  of  solute  (mg/liter)  are  denoted 
by  C  and  flow  rates  (liters/hour)  are  denoted  by  Q ,  each  with  subscripts  corresponding  to  the  specific 
tissue.  The  concentration  of  solute  in  the  air  is  denoted  by  Cc,  and  is  a  specified  quantity.  The  variables 
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TCE  concentration  in  adipose  tissue 


Figure  2:  Model  simulation:  Concentration  in  time  of  unbound  TCE  inside  the  perfusion-limited  fat  com¬ 
partment. 


Cvk,Cvm,Cvf,Cvbr,Cvi  and  Cvt  are  the  concentrations  of  solute  leaving  the  respective  organ  and  entering 
the  venous  blood  system.  In  this  case,  all  compartments  except  for  the  lung  are  perfusion-limited,  so  the 
concentration  of  solute  leaving  each  of  these  compartments  is  equal  to  the  concentration  of  free  solute  in 
that  compartment  itself.  In  the  kidney,  for  example,  this  implies 


where  Ck  is  the  total  concentration  of  solute  inside  the  kidney  compartment  and  P %  is  the  tissue/blood 
partition  coefficient  for  the  kidney. 

We  note  that  our  model  includes  a  venous  blood  compartment  governed  by  the  differential  equation  (2). 
In  this  type  of  PBPK  model,  the  algebraic  equation 

Cv  =  t y(Q  m^vm  T  Qt^vt  T  QfCvf  T  QbrCvbr  T  QlCyl  T  Qk^vk) 
h;c 

is  often  used  and  is  generally  sufficient  to  capture  the  behavior  of  the  solute  in  the  venous  blood.  This 
equation  is  based  on  a  steady-state  assumption  for  the  large  venous  blood  compartment.  Here  we  include 
the  more  complicated  model  (2)  because  it  will  be  more  compatible  with  our  modified  model  which  uses  a 
dispersion  compartment  for  the  adipose  tissue. 

Simulations  were  carried  out  using  the  model  (1)  -  (11)  to  predict  TCE  concentrations  in  Long- Evans 
rats.  The  Matlab  code  for  this  model  and  the  model  parameters  were  provided  by  Evans  et  al.  [22].  Figure  2 
depicts  an  example  simulation  for  the  following  input  function:  2000  ppm  TCE  in  the  surrounding  air  for 
one  hour,  followed  by  0  ppm  TCE  in  the  surrounding  air  for  14  hours.  The  concentration  of  unbound  TCE 
inside  the  perfusion-limited  fat  compartment  is  given  in  Figure  2.  Note  the  rapid,  exponential  decay  of 
concentrations  following  t  =  1  hour.  This  simulation  illustrates  that  this  type  of  compartmental  model  may 
not  be  able  to  capture  the  slow  accumulation  and  release  of  TCE  in  the  adipose  tissue. 
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2.2  A  PBPK  model  for  TCE  with  a  diffusion-limited  adipose  tissue  compart¬ 
ment 

A  second  type  of  compartment  used  in  PBPK  modeling  is  the  diffusion-limited,  or  membrane-limited  com¬ 
partment.  Unlike  the  perfusion-limited  model,  the  blood  flow  rate  to  the  tissue  is  rapid  compared  to  the 
transport  of  the  compound  across  the  cell  membranes  [38].  This  compartment  is  divided  into  two  subcom¬ 
partments:  the  intracellular  space  and  the  extracellular  space,  which  includes  the  vascular  blood  and  the 
interstitial  space.  The  concentration  of  compound  within  each  of  these  subcompartments  is  in  equilibrium, 
and  the  equations  governing  transport  of  the  compound  through  the  tissue  are 

=  Q{Cin  —  Cout)  +  —  Ce) 

=  KCs-Cj), 

where  Ve  and  Vj  are  the  volumes  of  the  extracellular  and  intracellular  compartments,  respectively,  and  Ce 
and  Ci  are  the  concentrations  of  compound  in  the  extracellular  and  intracellular  compartments,  respectively. 
Here  /x  is  the  cell  membrane  permeability  for  the  compound  (in  liters/hour),  and  Q,  Cin,  and  Cout  are  defined 
as  above. 

Now  we  modify  the  previous  PBPK  model  for  TCE  to  include  a  diffusion-limited  fat  compartment  (see 
Figure  3).  The  remainder  of  the  model  is  identical,  but  for  fat  we  assume  that  the  flow  of  blood  through 
adipose  tissue  is  much  faster  than  the  diffusion  of  solute  across  the  cell  membranes.  Thus,  we  replace  (1) 
by 

Vfe^lt  =  Qf(C*-Cc)+»(Cfi-Cfe)  (12) 

Vfi =  n(Cfe-Cfi)>  (13) 

where  Cfe  and  Cfi  are  the  concentrations  of  solute  in  the  extracellular  and  intracellular  spaces  of  the  fat 
tissue,  respectively,  Vfe  and  Vfi  are  the  volumes  of  these  extracellular  and  intracellular  spaces,  and  /x  is 
the  cell  membrane  permeability  coefficient  for  the  solute.  Our  new  system  of  equations  is  therefore  given 
by  (12),  (13)  and  (2)  -  (11). 

Model  simulations  were  carried  out  for  this  model  as  before.  The  concentrations  of  TCE  in  the  adipose 
subcompartments  are  shown  in  Figure  4  for  varying  values  of  /x.  For  /x  =  1 ,  the  concentration  of  unbound 
TCE  in  the  extracellular  subcompartment  as  seen  in  Figure  4  (Right)  is  nearly  identical  to  the  concentration 
of  unbound  TCE  in  the  perfusion-limited  adipose  compartment  (Figure  2). 

As  /x  decreases,  the  peak  concentration  of  unbound  TCE  in  the  intracellular  subcompartment  decreases 
and  the  decay  following  exposure  is  less  rapid.  The  concentration  profile  in  the  intracellular  space  for  /x  =  0.1 
shows  a  steady  increase  in  TCE  concentrations  before  t  =  1  hour,  followed  by  a  slow  decrease.  The  dynamics 
for  this  value  of  /x  appear  to  match  more  closely  with  the  expected  behavior  of  TCE  in  the  adipose  tissue 
than  the  dynamics  for  the  perfusion-limited  compartmental  model. 

One  disadvantage  with  the  diffusion-limited  model,  however,  is  that  there  is  only  one  extra  degree  of 
freedom  and  there  is  little  flexibility  in  obtaining  both  the  correct  concentration  profile  and  the  correct  peak 
value  in  the  intracellular  compartment.  Moreover,  since  TCE  has  a  low  molecular  weight  (131.3889)  and  is 
highly  lipophilic,  it  diffuses  easily  across  cell  membranes.  This  diffusion  rate  will  be  much  greater  than  the 
blood  flow  rate  to  the  adipose  tissue,  which  suggests  that  the  diffusion-limited  model  may  not  be  physically 
appropriate. 
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PBPK  Model  for  Inhaled  TCE 


Figure  3:  PBPK  model  for  inhaled  solute  with  membrane-limited  compartment  for  adipose  tissue. 
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Figure  4:  Model  simulation:  Concentrations  of  TCE  in  time  in  the  diffusion-limited  fat  compartment  for 
various  values  of  the  permeability  coefficient  fj,.  Left:  Total  concentration  of  TCE  in  the  extracellular 
subcompartment.  Right:  Concentration  of  unbound  TCE  in  the  intracellular  subcompartment. 


3  The  physiology  of  adipose  tissue  and  the  dispersion  model 


Neither  the  perfusion-limited  nor  the  diffusion-limited  compartmental  models  are  designed  to  account  for 
spatial  heterogeneities  within  the  intracellular  space  of  a  given  tissue.  Adipose  tissue  is  known  to  be  widely 
heterogeneous  in  both  its  composition  and  functional  properties.  Its  general  structure  includes  a  rich, 
irregular  vascular  bed,  with  capillaries  in  contact  with  each  adipocyte,  or  fat  cell  [47].  The  adipocytes, 
relatively  large  spherical  cells  containing  lipid  droplets  that  make  up  the  bulk  of  each  cell,  are  organized  into 
a  structure  of  lobules  supported  by  connective  tissue  [47].  See  Figure  5  for  an  illustration  of  a  representative 
adipocyte  and  capillary  unit,  which  is  suspended  in  the  surrounding  interstitial  fluid.  An  example  scanning 
electron  micrograph  of  adipose  tissue  that  is  on  the  same  scale  as  the  illustration  can  be  found  in  [49],  p.  92. 

The  major  metabolic  functions  of  adipose  tissue  include:  the  synthesis  and  storage  of  lipid  in  the  form 
of  triglyceride  [17],  the  breakdown  of  triglyceride  and  the  mobilization  of  lipid  in  the  form  of  free  fatty  acids 
[17],  and  the  uptake  of  fatty  acids,  glucose  and  amino  acids  [47]. 

Unlike  many  organs,  the  adipose  “organ”  is  located  in  multiple  distinct  depots  throughout  the  body. 
Several  studies  have  shown  that  there  are  pronounced  differences  in  tissue  composition,  blood  flow  and 
metabolic  activities  across  these  depots  [17],  [19],  [20],  [37],  [51],  which  suggest  that  a  single  perfusion-limited 
or  diffusion-limited  compartment  may  be  insufficient  to  capture  the  behavior  of  adipose  tissue  as  a  whole.  A 
possible  modification  to  these  models  is  the  inclusion  of  many  separate  perfusion-limited  or  diffusion-limited 
compartments,  each  representing  a  distinct  adipose  depot  and  each  with  its  own  physiological  parameters. 

In  addition  to  inter-depot  differences  in  blood  flow,  metabolism  and  tissue  composition,  there  are  hetero¬ 
geneities  within  each  adipose  depot  itself.  The  sizes  of  individual  fat  cells  greatly  vary,  as  does  the  amount 
of  lipid  in  each  adipocyte,  creating  an  uneven  distribution  of  lipids  across  the  tissue  [32].  Since  TCE  and 
other  compounds  of  interest  are  highly  lipophilic,  an  uneven  distribution  of  these  compounds  may  result, 
leading  to  a  spatially  dependent  concentration  function  inside  the  adipose  tissue. 

It  has  been  shown  that  the  metabolic  activities  of  adipose  tissue  are  directly  linked  to  its  blood  flow 
properties  [17].  As  the  metabolic  process  varies  in  different  regions  of  tissue  and  in  time,  the  perfusion 
of  blood  to  those  regions  changes  to  respond  to  the  needs  of  the  tissue  at  that  moment.  Moreover,  the 
blood  flow  to  adipose  tissue  is  directly  affected  by  the  local  concentrations  of  substrates  and  hormones  used 
in  the  process  of  metabolizing  lipids  [19].  Blood  flow  in  adipose  tissue  is  also  controlled  by  the  central 
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Figure  5:  Representation  of  an  adipocyte  and  capillary,  surrounded  by  interstitial  fluid. 


nervous  system,  which  can  alter  vascular  permeability,  capillary  density,  diffusion,  and  the  rate  of  transport 
of  lipophilic  compounds  [42].  Together,  these  physiological  characteristics  of  adipose  tissue  suggest  that 
a  spatially  varying  model  may  be  more  appropriate  for  describing  the  disposition  of  lipophilic  xenobiotics 
through  adipose  tissue. 


3.1  Modeling  philosophy 

An  important  feature  of  adipose  tissue  is  that  each  adipocyte  and  connecting  capillary  is  an  individual 
functional  “unit”  within  the  tissue.  There  are  millions  of  fat  cells  within  adipose  tissue,  and  it  would  be 
nearly  impossible  to  capture  the  disposition  of  TCE  specifically  through  each  and  every  adipocyte.  Rather, 
there  are  several  “aggregate”  modeling  approaches  that  may  be  used  to  approximate  the  transport  through 
the  tissue  without  actually  modeling  each  cell. 

One  such  approach  utilizes  mathematical  homogenization  theory,  as  used  in  [45]  and  [46]  to  estimate  the 
magnetic  permeability  of  magnetorheological  fluids.  The  homogenization  technique  involves  a  domain  with 
a  periodic  microstructure,  such  as  a  collection  of  individual  fat  cells  with  capillaries  immersed  in  interstitial 
fluid,  resulting  in  an  “effective”  equation  that  approximates  the  overall  behavior  inside  the  domain.  In  the 
case  of  adipose  tissue,  such  a  model  would  approximate  the  effective  concentration  of  TCE  inside  the  tissue. 

Another  modeling  approach  designed  to  handle  this  type  of  heterogeneity  makes  use  of  probability  theory. 
As  there  are  widely  varying  ranges  for  parameters  such  as  fat  cell  size,  blood  flow  rates  and  cell  permeability 
coefficients  inside  the  fat  tissue,  we  could  think  of  each  of  these  quantities  as  a  random  variable.  Each  of 
these  random  variables  would  have  an  unknown  probability  distribution  that  could  be  estimated  with  inverse 
problem  techniques.  See  [3]  and  [4]  for  detailed  descriptions  of  this  modeling  philosophy. 

A  third  type  of  model  sometimes  used  in  physiological  modeling  is  the  axial  dispersion  model,  which 
is  designed  to  account  for  the  heterogeneity  of  tissue.  It  assumes  that  the  amount  of  time  it  takes  for  a 
compound  to  pass  through  a  given  tissue  is  not  constant  for  each  particle  of  compound.  The  dispersion 
model  is  based  on  the  premise  that  the  variations  in  these  “residence  times”  are  a  result  of  radial  variations 
in  flow  velocity,  variations  in  path  length,  convective  mixing  in  the  direction  of  blood  flow,  and  molecular 
diffusion. 

One  important  feature  of  this  model  is  that  it  is  an  aggregate  type  of  model,  using  one  representative 
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“cell”  to  capture  the  behavior  of  a  large  collection  of  similar  “cells”  that  have  varying  properties.  The 
variability  among  these  cells  is  accounted  for  in  the  model  with  an  axial  dispersion  term,  and  the  axial 
dispersion  coefficient  is  a  representative  measure  of  the  variability  within  the  overall  system.  See  [35]  for  a 
detailed  development  of  the  axial  dispersion  model. 

Roberts  and  Rowland  [41]  adapted  this  type  of  dispersion  model  for  the  flow  of  compounds  through  the 
liver.  Based  on  this  model,  Banks,  Musante  and  Tran  [6]  developed  a  dispersion  model  for  the  transport 
and  elimination  of  2,3,7,8-tetrachlorodibenzo-p-dioxin  in  the  liver.  The  complex,  heterogeneous  physiology 
of  the  liver  and  the  widely  varying  circulation  of  blood  through  the  liver’s  sinusoids  both  suggest  this  type 
of  model,  and  simulations  carried  out  in  [5]  illustrate  the  effectiveness  of  this  model  for  describing  transport 
through  the  liver. 


4  A  dispersion  model  for  TCE  in  adipose  tissue 

Here  we  develop  an  axial  dispersion- type  model  for  the  adipose  tissue  compartment,  which  will  connect  with 
the  remaining  PBPK  compartmental  model.  This  model  is  based  on  the  dispersion  model  found  in  [35]  and 
will  be  adapted  for  a  representative  geometry  for  adipose  tissue. 


4.1  Geometry 

As  we  have  discussed,  adipose  tissue  is  composed  primarily  of  adipocytes,  which  are  relatively  large  spherical¬ 
shaped  cells.  The  adipocytes  are  supported  with  a  network  of  connective  tissue,  and  each  fat  cell  is  in  contact 
with  one  or  more  capillaries.  The  interstitial  space  makes  up  about  10%  of  the  total  tissue  weight  of  adipose 
tissue. 

Based  on  these  histological  characteristics,  we  choose  a  representative  geometry  for  the  fat  compartment 
that  includes  three  subcompartments:  the  capillary  subcompartment,  the  adipocyte  subcompartment  and 
the  interstitial  fluid  compartment.  We  choose  to  represent  the  adipocyte  subcompartment  by  a  single  sphere 
centered  at  the  origin  (see  Figure  6).  The  adipocyte  is  in  contact  with  a  capillary,  represented  by  a  cylindrical 
tube  with  a  circular  cross  section.  The  space  surrounding  the  capillary  and  the  adipocyte  is  occupied  by  the 
interstitial  fluid,  including  the  space  between  the  capillary  and  adipocyte. 

The  adipocyte  region  is  denoted  by  (A),  and  consists  of  the  sphere  centered  at  the  origin  with  radius  r±. 
The  capillary  or  blood  region  (B)  is  a  cylindrical  tube  with  a  circular  cross-section  that  wraps  around  the 
capillary.  This  cylindrical  tube  (again  see  Figure  6)  has  a  central  axis  r  =  r-2  +  ot,  9  =  9o,  £i  <  <j>  <  7r  —  £2  in 
spherical  coordinates,  where  r-2  >  7~i  >  0  and  E\,E2  >  0.  The  cross  section  of  the  capillary  has  radius  a  and 
area  Ab-  The  domain  of  this  cylindrical  tube  is  denoted  by  fls-  The  capillary  connects  with  the  arterial 
blood  at  <p  =  Ei,  and  similarly  connects  with  the  venous  blood  at  <j>  =  n  —  £2.  The  interstitial  fluid  region 
(I)  is  the  space  surrounding  the  adipocyte  and  the  capillary,  including  the  space  between  the  capillary  and 
adipocyte  (n  <  r  <  r?). 


4.2  Model  assumptions 

Blood  flows  from  the  arterial  compartment  into  the  capillary,  carrying  TCE  to  the  adipose  tissue.  TCE  then 
diffuses  back  and  forth  across  the  capillary  endothelium  into  and  out  of  the  interstitial  fluid,  and  TCE  in 
the  interstitial  fluid  may  also  diffuse  back  and  forth  across  the  plasma  membrane  of  the  adipocyte  into  and 
out  of  the  fat  cell.  At  the  exit  to  the  capillary  the  blood  and  TCE  flow  into  the  venous  blood  compartment. 
See  Figure  6  for  a  schematic  of  the  assumed  geometry. 
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Figure  6:  Adipose  tissue  represented  geometrically:  the  adipocyte  region  (A)  is  a  sphere,  surrounded  by  the 
interstitial  space  (I).  The  capillary  or  blood  region  (B)  is  a  cylindrical  tube  that  wraps  around  the  capillary. 
Coordinates  are  in  spherical  coordinates  (r,9,<j>). 
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Some  studies  have  suggested  that  there  may  be  an  “interfacial  continuum”  between  capillary  endothelial 
cells  and  adipocytes  [9],  [10],  [11],  [44].  That  is,  there  is  evidence  of  membrane  continuity  between  the 
membranes  of  capillary  endothelial  cells  and  the  plasma  membranes  of  adipocytes.  Since  TCE  is  highly 
lipophilic  and  the  cell  membranes  are  composed  of  lipid  bilayers  [31],  TCE  may  diffuse  directly  across  the 
capillary  endothelium  and  the  adipocyte  plasma  membrane  at  these  areas  of  membrane  continuity.  This 
would  result  in  direct  transport  between  the  capillary  and  the  adipocyte. 

Therefore  there  are  two  feasible  routes  of  transport  between  the  capillary  and  the  adipocyte:  indirect 
transport  via  diffusion  through  the  interstitial  fluid  and  direct  transport  across  an  “interfacial  continuum.” 
It  would  be  difficult  to  determine  in  vivo  the  exact  proportions  of  direct  and  indirect  transport  of  TCE 
occurring  between  an  adipocyte  and  a  capillary.  Our  model  will  include  both  of  these  types  of  transport, 
and  the  proportions  of  each  type  will  be  parameters  that  may  be  estimated  using  inverse  problem  techniques. 

We  note  that  here  we  do  not  attempt  to  include  in  our  representative  geometry  the  areas  of  continuity 
between  the  capillary  endothelium  and  the  adipocyte  plasma  membrane.  Rather  than  modeling  each  of  these 
regions  individually,  we  will  include  a  general  term  in  our  model  that  will  account  for  any  direct  transport 
that  occurs. 

We  assume  that  the  concentrations  of  TCE  in  the  interstitial  fluid  are  in  equilibrium  at  points  far  away 
from  the  adipocyte  and  the  capillary.  Specifically,  we  assume  that  the  net  radial  flux  of  TCE  across  the 
spherical  shell  r  =  r3  is  equal  to  zero  for  some  r3  >  r2  +2a,  where  r2  +2a  is  the  outer  radius  of  the  capillary. 
We  define  this  region  of  the  interstitial  space  as 

fij  =  {(r,  0.  (f>)  :  r'i  <  r  <  r3, 0  <  9  <  2n,  0  <  <j>  <  n}. 


Similarly,  since  the  adipocyte  is  large  relative  to  the  size  of  the  capillary  and  nearby  interstitial  fluid,  we 
assume  that  the  net  radial  flux  of  TCE  inside  the  adipocyte  across  the  spherical  shell  r  =  r0  is  equal  to  zero 
for  some  0  <  r0  <  r\.  We  thus  define  the  adipocyte  domain  by 

ff A  =  { 6,<f>)  :  r0  <  r  <  rq,  0  <  9  <  2 n,  0  <  <j>  <  n}. 

Here  we  assume  that  the  difference  between  the  inner  radius  r0  and  the  outer  radius  r3  are  small  compared 
to  the  radius  r\  of  the  adipocyte. 

The  remaining  assumptions  for  our  dispersion  model  are  based  on  those  used  by  Roberts  and  Rowland 
for  the  liver,  adapted  here  for  the  physiology  of  adipose  tissue. 


1.  Conditions  in  the  capillary  do  not  vary  in  a  cross  section  normal  to  flow.  This  implies  that  the 
concentration  Cb  of  TCE  in  the  capillary  depends  only  on  <j>.  The  basis  for  this  assumption  is  that 
the  flow  in  the  r  and  0  directions  is  negligible  compared  to  the  flow  of  the  blood  in  the  <p  direction. 


2. 


Axial  diffusion  of  solute  is  slow  compared  to  the  convection  of  solute  in  the  capillary,  so  it  may  be 
ignored.  It  can  be  shown  that  without  loss  of  generality,  a  separate  axial  diffusion  term  can  be  omitted 
in  any  axial  dispersion  model.  If  we  assume  axial  diffusion,  then  we  obtain  the  following  diffusive  mass 
flux  in  the  <p  direction: 


where  C  is  the  concentration  of  solute  and  D  is  the  axial  diffusivity  coefficient.  When  considering  axial 
dispersion,  the  dispersive  mass  flux  in  the  <j>  direction  is  given  by 


Ja  —  —V 


dC_ 

do'' 


where  V  is  the  axial  dispersivity  coefficient.  Since  neither  of  these  coefficients  are  known  in  general,  they 
must  be  estimated  using  inverse  problem  techniques.  Indeed,  the  dispersivity  coefficient  is  not  physical, 
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since  it  is  an  effective  coefficient  which  accounts  for  variations  in  residence  times,  path  length,  local 
geometry,  etc.  Therefore  we  can  interpret  the  coefficient  of  this  flux  term  as  either  the  axial  diffusivity, 
axial  dispersivity,  or  as  a  net  axial  transport  coefficient  involving  a  combination  of  the  two. 

3.  The  transport  of  TCE  between  the  capillary  and  the  interstitial  space  is  assumed  to  obey  Fick’s  law 
of  diffusion,  as  are  the  transport  between  the  interstitial  space  and  the  adipocyte  and  the  transport 
between  the  capillary  and  the  adipocyte. 

4.  No  metabolism  or  elimination  of  solute  occurs  in  adipose  tissue.  This  assumption  is  based  on  the  fact 
that  TCE  is  metabolized  by  the  enzyme  Cytochrome  P450  2E1,  which  is  located  primarily  in  the  liver 
and  is  not  known  to  be  present  at  significant  levels  inside  adipose  tissue.  Moreover,  TCE  is  eliminated 
only  via  exhalation  from  the  lungs  and  via  the  excretion  of  its  metabolites  in  the  urine. 

5.  Only  unbound  solute  can  cross  cell  membranes.  This  is  a  standard  assumption  used  in  physiological 
modeling. 

6.  Concentrations  of  TCE  in  the  interstitial  fluid  immediately  surrounding  the  capillary  and  adipocyte 
(ri  <  r  <  rff)  are  uniform  with  respect  to  r,  depending  only  on  9  and  <f>.  This  assumption  is  based  on 
the  small  radial  distance  between  r±  and  r 3. 

7.  Concentrations  of  TCE  in  the  interstitial  fluid  near  the  capillary  vary  only  slightly  with  respect  to  9, 
and  the  transport  between  the  capillary  and  the  interstitial  space  is  approximated  by  the  transport 
occurring  at  r  =  r2,  9  =  9o,  £1  <  <j>  <  n  —  e2. 

8.  Concentrations  of  TCE  inside  the  adipocyte  near  the  outer  radius  r±  are  uniform  with  respect  to  r, 
depending  only  on  9  and  <f>. 

9.  The  ratio  fs  of  the  unbound  solute  to  the  total  solute  concentration  within  the  blood  is  a  constant.  This 
is  also  a  standard  physiological  assumption,  based  on  the  fact  that  the  total  and  unbound  concentrations 
of  solute  in  the  blood  are  in  equilibrium  with  each  other,  and  a  linear  relationship  between  the  free 
and  unbound  concentrations  is  assumed.  Let  Cb  and  CUB  denote  the  blood  concentrations  of  total 
and  unbound  solute,  respectively.  Then 

Cub  =  JbCb 

for  some  /#,  and  the  fact  that  these  concentrations  are  in  equilibrium  implies  that  fs  must  be  a 
constant. 

10.  The  ratio  //  of  the  unbound  solute  to  the  total  solute  concentration  within  the  interstitial  space  is  a 
constant.  This  assumption  is  the  interstitial  space  analog  of  Assumption  9. 

11.  The  ratio  Ja  of  the  unbound  solute  to  the  total  solute  concentration  within  the  adipocytes  is  a  constant. 
This  assumption  is  the  adipocyte  analog  of  Assumption  9. 


4.3  Derivation  of  equations:  Capillary  region 


Consider  a  “spherical  cuboid”  volume  element  spanning  from  r  to  r  +  Ar,  9  to  9  +  A 9  and  (j>  to  </>  +  A <f>  (see 
Figure  7).  The  volume  of  this  element  is  Ay  =  r2  sin(j)ArA(j)A9,  and  we  will  denote  its  surface  by  AS. 

In  the  capillary,  the  concentration  of  TCE  is  uniform  with  respect  to  both  r  and  9,  and  we  have  axial 
dispersion  in  the  <j>  direction.  Moreover,  we  assume  that  a  is  small  relative  to  r2  so  that  we  may  approximate 
r  w  r2  for  r  £  fl#.  This  implies  that  the  flux  vector  J  =  ( Jr,Jg ,  Jq)  is  given  by 


J  =  (  0,0, 


'Db  9Cb 
r2  dcj> 


+  vCb), 
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Figure  7:  “Spherical  cuboid”  volume  element. 


where  T>b  is  the  axial  dispersion  coefficient  with  units  m2/hr.  The  velocity  v  is  given  by 

_  Qf  —  r  ( q  litersN\  (  1  f  1  meter3  \ 

V  1000AB  V  ^  hour )  meters2/  \  1000  liters/ 

with  units  m/hr,  where  <5/  is  the  blood  flow  rate  to  the  fat  tissue  in  liters/hr  and  Ab  is  the  cross-sectional 
area  of  the  capillary  in  meters2 . 

The  dimensionless  coefficient  T  is  included  to  account  for  the  inversion  of  the  Fahraeus-Lindquist  effect, 
which  is  a  phenomenon  that  involves  the  flow  of  blood  as  a  function  of  the  vessel  diameter.  To  understand 
the  phenomenon  involved,  consider  the  flow  of  blood  from  a  larger  vessel  to  a  smaller  vessel,  with  the  smaller 
vessel  ranging  between  15  to  500  microns  in  diameter.  If  we  assume  that  blood  with  constant  hematocrit  (the 
fraction  of  blood  volume  made  up  of  the  red  blood  cells)  flows  from  the  larger  vessel  to  the  smaller  vessel,  the 
hematocrit  in  the  smaller  vessel  decreases  as  the  diameter  of  the  vessel  decreases  [27].  This  is  known  as  the 
Fahraeus  effect,  and  is  related  directly  to  the  Fahraeus-Lindquist  effect,  which  involves  a  concurrent  decrease 
in  the  apparent  viscosity  of  blood  in  the  smaller  vessel  as  the  diameter  and  the  hematocrit  decrease  [27]. 

The  apparent  viscosity  fj,app  can  be  expressed  as  a  function  of  the  vessel  radius  R  and  length  L.  the  drop 
in  pressure  A P  across  the  length  of  the  vessel,  and  the  volumetric  flow  rate  Q  as  in  [26]  and  [27]: 

_  nR4AP 
^app  ~  8 LQ  ' 

A  different  phenomenon  occurs  in  vessels  less  than  about  4-6  microns  in  diameter.  These  vessels  are 
significantly  smaller  than  the  normal  8  micron  diameter  of  red  blood  cells,  making  it  difficult  for  the  red 
blood  cells  to  enter  vessels  of  this  size  [27].  This  results  in  an  increase  of  the  flow  resistance  and  hence  an 
increase  in  the  apparent  viscosity  of  the  blood  as  the  diameter  of  the  vessel  decreases,  which  is  the  inversion 
of  the  Fahraeus-Lindquist  effect  [27].  This  phenomenon  leads  to  a  concurrent  decrease  in  the  flow  rate  of 
the  blood  inside  the  smaller  tube  compared  to  the  flow  rate  in  the  larger  feed  vessel. 

The  capillaries  in  adipose  tissue  typically  have  lumen  diameters  measuring  between  2  and  4  microns  [19], 
which  would  place  them  in  the  range  of  vessels  covered  by  the  inversion  of  the  Fahraeus-Lindquist  effect.  This 
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implies  a  significant  difference  between  the  flow  velocities  of  blood  inside  the  arterial  blood  system  versus 
the  velocities  inside  the  capillaries  of  the  adipose  tissue.  Hence  we  include  the  “reverse  Fahraeus-Lindquist 
parameter”  T  <  1  to  account  for  this  decrease  in  blood  flow  rates  as  the  blood  and  TCE  flow  from  the 
arterial  vessels  to  the  adipose  capillaries. 

Conservation  of  mass  implies  that 


Net  rate  of  change 
of  TCE  in  blood  element 


+ 


(  Net  flow  rate  of  TCE 
within  blood  element 


/  Net  flow  rate  of  TCE 
[  between  blood  element 
\  and  adipocyte  region 


Net  flow  rate  of  TCE 
between  blood  element 
and  interstitial  space 


First  we  consider  only  transport  from  within  the  capillary  (i.e.,  no  exchange  with  the  interstitial  space 
or  adipocyte  region).  Then 

(  Net  flow  rate  of  TCE  \  _  f  f  „ 

V  within  blood  element  J  J  JAS  n  ’ 

where  n  is  the  outer  unit  normal  to  the  surface  AS.  Dividing  by  the  volume  AV  of  the  element,  we  obtain 

=  A lvIL-’'”dS 
=  Zrffl^-W 


Net  flow  rate  of  TCE 
within  blood  element 
per  unit  volume 


by  the  Divergence  Theorem.  Taking  the  limit  as  AV  — >  0  at  r  =  we  have 


(  Net  flow  rate  of  TCE  \ 
in  the  blood  element 
per  unit  volume 
\  at  the  point  (r2,6,4>)  J 


Div(-J) 

-iiir'i Ur)~ 


1  djg 


1 


r,  dr 

1 


r2  sin  <p  89  r 2  sin  <fi  1 


-(sin  (f)J<p) 


8 


r2  sin  (j)  84 > 

1  8 


r2  sin  4>  84> 

Multiply  both  sides  by  the  volume  of  the  capillary  Vb  to  obtain 


(sin  4 Wf) 

VBdCB  ,  „ 

sin  cp  ( - +  vC b 


r2  84> 


Net  flow  rate  of  TCE 
within  capillary 


Vb  8 
r 2  sin  4>  84> 


.  ,  VB  8Cb  r, 

sin  4> - — vCB 

r2  84> 


Now  consider  the  exchange  of  TCE  between  the  capillary  and  interstitial  space.  By  Assumption  3,  this 
transport  is  governed  by  Fick’s  law  of  diffusion,  so  that 


^  =  XlHBl(flCl(do )  -  f bCb ), 


Net  flow  rate  of  TCE 
between  blood  element 
and  interstitial  space 
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where  /ibi  is  the  permeability  coefficient  for  transport  between  the  capillary  and  interstitial  regions,  with 
units  liters/hour.  The  dimensionless  coefficient  A j  represents  the  fraction  of  transport  that  occurs  between 
the  capillary  and  interstitial  fluid  (versus  the  direct  transport  between  the  capillary  and  the  adipocyte). 


Similarly,  we  model  the  transport  between  the  capillary  and  the  adipocyte  by 
/  Net  flow  rate  of  TCE  \ 

between  blood  element  =  ^a^baUaCa(9q)  ~  IbCb), 
\  and  adipocyte  J 


where  Hba  is  the  permeability  coefficient  for  transport  between  the  capillary  and  adipocyte  regions  (in 
liters/hour),  and  A  .4  represents  the  fraction  of  transport  that  occurs  between  the  capillary  and  the  adipocyte. 
Note  that  we  have  Aj  +  A^  =1. 


Since  the  total  rate  of  change  of  TCE  at  a  point  in  the  blood  element  is  given  by  we  have 


VB 


9Cb  _  Vb  d 
dt  r 2  sin  cf>  d(f> 


sin< 


Vb  9Cb 
r2  dcj) 


■vCB 


+  Ai/iB/(/iCj(#o)  -  IbCb)  +  ^aHbaUaCa(Oq)  -  f bCb ) 


for  t  >  0  and  (r,  6,<j>)  £  fi#,  where 

&b  =  {( r,9 ,  <j>)  :  r  =  r2,9  =  60,£i  <  <t>  <  n  -  £2}  C  fis- 


Note  that  Cb  =  Csit,  4>)  due  to  uniformity  in  r  and  9. 


4.3.1  Boundary  conditions 


To  obtain  boundary  conditions  for  the  capillary  region,  we  use  flux  balance  principles.  The  mass  flux  in 
the  capillary  across  the  boundary  (f>  =  £1  is  given  by  J<p\rp=£1 ,  and  the  mass  flux  in  the  arterial  blood  across 
<f>  =  £1  is  (<2c/(1000Ab))Co.  Therefore  we  have 


Vb  QCb 
r2  dip 


(t,  <j>)  +  vCB(t,  <j>) 


<t>=£  1 


Qc 

1000 AB 


Ca(t). 


Similarly,  using  flux  balance  at  the  boundary  (j>  =  ir  —  s2  between  the  capillary  and  the  venous  blood,  we 
obtain 


Vb  9Cb 
r2  dcj) 


(' t ,  <t>)  +  vCB(t,  <t>) 


<j>= 7T  — £2 


Qc 

1000.4^ 


Cv(t). 


4.4  Interstitial  space 


In  the  interstitial  space,  TCE  enters  and  exits  via  diffusion  from  the  capillary  region  and  via  diffusion  from 
the  adipocyte. 


We  assume  that  the  interstitial  space  immediately  surrounding  the  adipocyte  and  the  capillary  has  a 
radial  thickness  sufficiently  small  so  that  diffusion  in  the  radial  direction  is  nearly  instantaneous,  i.e.,  the 
concentration  of  TCE  in  the  interstitial  space  depends  on  (f ,  9,  <f>)  but  not  r.  Therefore,  we  may  assume  that 
r  m  r\  in  the  interstitial  space.  Inside  the  interstitial  fluid,  TCE  diffuses  across  the  interstitial  space  in  both 
the  9  and  <j>  directions. 


Based  on  these  assumptions,  the  mass  flux  vector  for  the  interstitial  space  is  given  by 


J  —  (Jr,  Jfj ,  Jrj,) 


PT  d_Cj_ 
ri  sin  <j)  99 


Vi  9CA 
n  9(f)  )  ' 
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The  parameter  £>/  is  the  standard  diffusion  coefficient  with  units  m2/hr. 

Now  consider  an  arbitrary  “spherical  cuboid”  volume  element  in  the  interstitial  space,  as  depicted  in 
Figure  7.  Using  mass  balance  as  before,  we  have 


Net  rate  of  change  of  TCE 
in  interstitial  element 


+ 


f  Net  flow  rate  of  TCE  \ 
V  within  interstitial  element  J 


/  Net  flow  rate  of  TCE  \ 

I  between  interstitial  + 
y  element  and  capillary  J 


Net  flow  rate  of  TCE 
between  interstitial 
element  and  adipocyte 


First  we  consider  the  transport  of  TCE  within  the  interstitial  element  only  (i.e.,  no  exchange  with  the 
other  regions).  Then 

Net  flow  rate  of  TCE  \ 
within  interst.  element  J 

where  n  is  the  outer  unit  normal  to  the  surface  AS.  Dividing  by  the  volume  AU  of  the  element,  we  obtain 


- 

J  Jas 


—  J  ■  n  dS , 


/  Net  flow  rate  of  TCE 
I  within  interst.  element 
\  per  unit  volume 


—  J  •  n  dS. 


Using  the  Divergence  Theorem  and  taking  the  limit  as  AU  — >  0  at  r  =  n,  we  have 


(  Net  flow  rate  of  TCE  \ 
within  interst.  element 
per  unit  volume 
\  at  the  point  (n  ,8,<p)  J 
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sin  <j>- 


d(p 


As  with  the  capillary,  we  multiply  by  the  volume  and  then  add  in  the  transport  between  the  other  regions 
to  obtain 


,  9Cj 
1  0 1 


VjDi  r  i  o-Ci  i  o  ( .  QCiX 
r2  _sin2^  082  +  sirup  d(f>  \  dcp  J 
+  &e0 {9)xb{<P)Xhbi{SbCb  ~  fiCj)  +  HiaUaCa  ~  fiCi) 


(14) 


for  t  >  0  and  (r,  9,  cp)  €  fij,  where 


&i  =  {(r,9,(p)  :  r  =  r\ ,  0  <  9  <  27t,  0  <  <p  <  n}  C  D/ 


and  xb  is  the  following  indicator  function  in  cp: 

xBU)  =  {  1  if  <^e(£i^-£2) 

A  w  (0  otherwise, 

This  function  xb  is  used  with  dg0,  the  Dirac  delta  distribution  centered  at  9  =  9o,  to  account  for  the  capillary 
-  interstitial  space  transport  (and  the  capillary  -  adipocyte  transport) . 
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4.4.1  Boundary  conditions 


Since  we  have  diffusion  in  the  interstitial  space  in  both  the  0  and  9  directions,  we  require  boundary  conditions 
in  0  and  9.  The  geometry  of  the  spherical  shell  dictates  periodic  boundary  conditions  in  9.  That  is, 

t  >  0  and  0  <  0  <  n 

t  >  0  and  0  <  0  <  n. 

If  we  carry  out  separation  of  variables  on  the 
PDE  (14)  for  the  interstitial  space,  we  obtain  an  equation  in  0  which  is  a  type  of  singular  Sturm-Liouville 
problem  (see  [16],  p.  324).  These  problems  have  the  form 

( pu')'  —  qu  +  \pu  =  0, 

where  the  variable  coefficient  p  vanishes  at  one  or  both  of  the  endpoints.  In  our  case,  we  have  p(<j>)  =  sin  0, 
which  vanishes  at  both  0  =  0  and  0  =  7 r. 

The  corresponding  boundary  value  problems  require  a  finiteness  or  continuity  condition  at  the  singular 
endpoints,  taking  the  place  of  the  usual  type  of  boundary  condition  [16].  Here  we  will  require  the  finiteness 
conditions 


Ci(t,9,4 f>)  =  Ci(t,9  +  2n,  0), 
^M,^)  =  ^M  +  2tt,0), 


The  boundary  conditions  in  0  are  less  straightforward. 


Ci(t,9,0)  <  oo 

Ci(t,9, 7r)  <  oo. 

4.5  Derivation  of  equations:  Adipocyte  region 

In  the  adipocyte  region,  TCE  enters  and  leaves  the  fat  cell  via  diffusion  across  the  cell  membrane  between  the 
adipocyte  and  the  interstitial  space,  and  via  diffusion  between  the  adipocyte  and  capillary.  We  assume  that 
diffusion  of  TCE  in  the  radial  direction  inside  the  adipocyte  is  nearly  instantaneous,  so  that  the  concentrations 
are  uniform  in  r  for  ro  <r  <r\.  As  with  the  interstitial  space,  TCE  in  the  adipocyte  diffuses  across  the  fat 
cell  in  the  9  and  0  directions. 

Using  flux  balance  principles  as  before,  we  obtain  the  equation  (taking  the  limit  AV  — »  0  at  r  =  r0) 

T  f  dCA  _  VaDa 
A  &t  r2 

+  $00  {6)XB{4>)^AHBA{f bCb  -  SaCa)  +  fiiAifiCj  -  f aC a)  (15) 

for  t  >  0  and  (r,  9, 0)  £  &a,  where 

&A  =  {(r,  9,  0)  :  r  =  r0, 0  <  9  <  0  <  0  <  n}  C 

As  before,  the  parameter  DA  is  the  diffusion  coefficient  with  units  m2/hr. 


1  d^CA  1  d  f  .  dCA\ 
sin2  0  d9 2  sin  0  90  V  90  ) 


4.5.1  Boundary  conditions 

The  boundary  conditions  for  the  adipocyte  region  are  similar  to  those  for  the  interstitial  space,  and  are 
derived  using  the  same  principles.  We  require  periodic  boundary  conditions  in  9  and  finiteness  conditions  in 
0: 

t  >  0  and  0  <  0  <  7r 
t  >  0  and  0  <  0  <  7r, 


CA(t,  9, 0)  =  CA(t,9  +  2tt,<P), 
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Ca  {t,  8,0)  <  oo 

CA{t,0,n)  <  oo. 


4.6  Summary  of  model  equations 

Combining  the  PBPK  whole-body  model  with  the  adipose  submodel,  we  obtain 
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Vi  — 

1  dt 

=  Ql(Ca  Cvl) 

+  C' vi 

(35) 
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dt 
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4.7  Model  simulations 

We  used  the  TCE  adipose  dispersion  model  (16)  -  (28)  coupled  with  the  whole-body  PBPK  model  (29)  -  (38) 
to  generate  model  simulations  with  the  same  input  function  and  non-adipose  parameters  as  in  Sections  2.1 
and  2.2.  In  these  simulations  we  have  chosen  to  approximate  ro  =  r±  =  The  basis  for  this  approximation 
is  the  assumption  that  the  radial  distance  between  ro  and  r$  is  small  compared  to  the  radius  r±  of  the 
adipocyte.  Moreover,  the  radial  variable  r  occurs  in  our  system  of  equations  only  as  a  coefficient,  which 
is  fixed  at  r  =  r2  in  (16),  r  =  n  in  (19),  and  r  =  ro  in  (24).  Therefore  the  approximation  ro  =  r\  =  r2 
in  (16)  -  (28)  should  have  little  impact  on  the  resulting  solutions. 

Figure  8  depicts  an  example  simulation  for  the  concentration  of  unbound  TCE  inside  the  adipocyte  region 
versus  time  for  several  points  along  the  “equator”  <f>  =  7t/2,  where  we  have  chosen  9 0  =  n.  Note  that  these 
concentration  profiles  look  similar  to  those  for  the  diffusion-limited  model  in  Figure  4. 

An  important  difference  between  the  diffusion-limited  model  simulations  and  the  dispersion  model  simu¬ 
lations  is  that  all  of  the  concentration  profiles  in  Figure  8  were  generated  with  fixed  permeability  parameters 
(Hba  =  0.03,  ubi  =  0.02  and  hia  =  0.025  liters/hour,  Aj  =  Xa  =  0.5)  and  observations  taken  at  varying 
locations,  while  those  shown  in  Figure  4  were  generated  with  varying  values  of  the  permeability  coefficient 
fi.  By  design,  the  dispersion  model  accounts  for  spatial  variation  in  concentrations  within  each  region  of  the 
adipose  tissue,  and  can  be  used  to  predict  various  concentration  profiles  within  the  fat. 

Moreover,  by  varying  the  adipose  model  parameters  T  and  Ab,  the  reverse  Fahraeus-Lindquist  parameter 
and  the  cross-sectional  area  of  the  capillary  respectively,  we  can  control  the  peak  concentration  values  inside 
the  adipose  tissue  regions.  This  is  accomplished  by  modifying  the  flux  between  the  adipose  capillary  and 
the  arterial  blood  compartment,  as  well  as  the  flux  between  the  adipose  capillary  and  the  venous  blood 
compartment.  The  concentration  profiles  in  Figure  8  were  generated  with  T  =  0.16  and  Ab  =  10  4,  and 
Figure  9  depicts  concentrations  for  T  =  0.08  and  Ab  =  5  x  10~6.  Note  that  for  a  given  value  of  6 ,  the 
concentration  curves  are  similar  in  shape  but  vary  in  their  peak  values. 

Concentration  profiles  also  can  be  altered  significantly  by  changing  the  proportions  Xa  and  A j  of  direct 
and  indirect  transport  between  the  capillary  and  the  adipocyte.  Figures  10  and  11  depict  simulations  of 
unbound  TCE  adipocyte  concentrations  along  < j>  =  7r/2  for  varying  values  of  Xa,  Aj.  In  Figure  10,  the  indirect 
transport  is  predominant  over  the  direct  transport  (A^  =  0,  Aj  =  1  and  A  a  =  0.4,  A  j  =  0.6),  while  the  direct 
transport  predominates  in  Figure  11  (Xa  =  0.6,  A  j  =  0.4  and  Xa  =  0.65,  A  j  =  0.35).  As  Xa  increases  and 
Aj  decreases,  the  levels  of  TCE  increase  and  the  shapes  of  the  curves  change  dramatically. 

As  seen  in  these  example  simulations,  the  dispersion  model  is  able  to  predict  concentration  curves  of 
TCE  inside  the  adipose  tissue  with  varying  levels  of  accumulation  and  rates  of  decay.  This  suggests  that  this 
model  might  be  successfully  used  with  parameter  estimation  techniques  to  predict  experimentally  determined 
dynamics  of  TCE  in  adipose  tissue.  A  paper  discussing  the  computational  efforts  introduced  here  and  the 
associated  parameter  estimation  problem  is  currently  under  preparation. 
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TCE  adipocyte  concentrations  along  4>  =  it/2 


Figure  8:  Model  simulations:  Concentrations  in  time  of  unbound  TCE  inside  the  adipocyte  region  of  the 
adipose  tissue  at  (f>  =  7t/2  and  9=0,  n/A,  7t/2,  37t/4  and  n,  with  J7  =  0.16,  Ab  =  10  4,  A j  =  0.5  and 
Xa  =  0.5. 


TCE  adipocyte  concentrations  along  §  =  nl2 


Figure  9:  Model  simulations:  Concentrations  in  time  of  unbound  TCE  inside  the  adipocyte  region  of  the 
adipose  tissue  at  <j>  =  7t/2  and  9  =  0,  n/A,  n/2,  37t/4  and  n,  with  J7  =  0.08,  Ab  =  5  x  10  5,  A j  =  0.5  and 
Xa  =  0.5. 
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TCE  adipocyte  concentrations  along  $  =  it/2 


TCE  adipocyte  concentrations  along  $  =  n/2 


Figure  10:  Model  simulations:  Concentrations  in  time  of  unbound  TCE  inside  the  adipocyte  region  of  the 
adipose  tissue  at  <f>  =  n/2  and  9  =  0,  7t/4,  7t/2,  3n/4  and  n,  with  T  =  0.16,  Ab  =  10  4,  and  Left:  A /  =  1, 
A  a  =  0;  Right:  A  /  =  0.6,  A  a  =  0.4. 


TCE  adipocyte  concentrations  along  <|>  =  ir/2  TCE  adipocyte  concentrations  along  <j>  =  jt/2 


Figure  11:  Model  simulations:  Concentrations  in  time  of  unbound  TCE  inside  the  adipocyte  region  of  the 
adipose  tissue  at  (j>  =  7t/2  and  6  =  0,  7t/4,  7t/2,  37t/4  and  n,  with  T  =  0.16,  Ab  =  10  4,  and  Left:  A /  =  0.4, 
A^  =  0.6;  Right:  A  /  =  0.35,  A  a  =  0.65. 
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5  Concluding  remarks 


In  this  paper  we  have  discussed  various  modeling  techniques  for  capturing  the  transport  behavior  of  trichloro¬ 
ethylene  in  the  body,  and  specifically  inside  the  adipose  tissue.  TCE  is  a  member  of  a  class  of  chemicals  that 
is  known  to  accumulate  inside  the  adipose  tissue,  and  its  persistence  there  has  a  major  impact  on  the  overall 
kinetics  of  TCE  inside  the  body.  Physiologically  based  pharmacokinetic  modeling  is  the  standard  approach 
for  describing  the  disposition  of  compounds  such  as  TCE  in  the  body,  and  is  based  on  a  system  of  coupled 
compartmental  equations  that  each  model  a  specific  tissue  or  organ. 

First  we  implemented  the  standard  PBPK  model  to  simulate  inhaled  TCE  in  Long-Evans  rats,  using 
a  perfusion-limited  compartment  for  the  adipose  tissue.  We  then  adapted  this  model  by  changing  the  fat 
compartment  to  a  diffusion-limited  model.  The  perfusion-limited  model  assumes  that  the  entire  tissue  is  in 
rapid  equilibrium  and  the  diffusion-limited  model  assumes  that  the  intracellular  space  and  extracellular  space 
of  a  tissue  are  individually  in  rapid  equilibrium.  These  “well-stirred”  assumptions  may  not  be  adequate  to 
capture  the  behavior  of  TCE  inside  the  fat  tissue. 

Adipose  is  a  significantly  heterogeneous  tissue,  with  widely  varying  physiological  properties  such  as  fat  cell 
size,  lipid  distribution,  blood  flow  rates  and  cell  permeabilities.  Each  adipocyte  is  a  relatively  independent 
functional  unit  with  its  own  blood  supply,  suggesting  that  a  lipophilic  compound  such  as  TCE  may  have 
more  complicated  transport  dynamics  in  the  fat  tissue.  These  dynamics  are  likely  to  be  spatially  varying, 
reflecting  the  spatial  heterogeneities  in  the  adipose  tissue  itself. 

To  account  for  these  heterogeneities  and  more  closely  model  the  specific  physiology  of  adipose  tissue,  we 
developed  an  axial  dispersion  type  model  to  describe  the  transport  of  TCE  inside  the  fat.  This  aggregate 
model  uses  a  representative  adipocyte-capillary  structure  to  approximate  the  dynamics  occurring  in  individ¬ 
ual  adipocytes  throughout  the  tissue.  The  representative  geometry  and  the  model  itself  are  based  on  the 
physiological  features  of  adipose  tissue  and  the  expected  behavior  of  lipophilic  compounds  in  the  body. 

Model  simulations  were  generated  for  various  parameter  sets,  demonstrating  that  this  type  of  model 
might  be  well-suited  to  predict  the  experimental  behavior  of  TCE  inside  adipose  tissue  using  parameter 
estimation  techniques.  Moreover,  these  simulations  demonstrate  that  concentration  profiles  depend  greatly 
on  the  choice  of  model  parameters  and  on  the  predominant  type  of  transport  between  the  capillary  and 
adipocyte  regions.  Based  on  these  results  and  comparisons  with  the  two  standard  models,  it  appears  that 
accounting  for  the  specific  physiology  of  adipose  tissue  is  important  in  developing  an  accurate  model  for  the 
transport  of  TCE  in  the  adipose  tissue  and  the  body. 
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